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Dynamics of influence processes on networks: 
Complete mean-field theory; the roles of response functions, connectivity, and 
synchrony; and applications to social contagion 
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Vermont Complex Systems Center, and Computational Story Lab, 
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We study binary state dynamics on a network where each node acts in response to the average 
state of its neighborhood. Allowing varying amounts of stochasticity in both the network and node 
responses, we find different outcomes in random and deterministic versions of the model. In the limit 
of a large, dense network, however, we show that these dynamics coincide. We construct a general 
mean field theory for random networks and show this predicts that the dynamics on the network are 
a smoothed version of the average response function dynamics. Thus, the behavior of the system can 
range from steady state to chaotic depending on the response functions, network connectivity, and 
update synchronicity. As a specific example, we model the competing tendencies of imitation and 
non-conformity by incorporating an off-threshold into standard threshold models of social contagion. 
In this way we attempt to capture important aspects of fashions and societal trends. We compare 
our theory to extensive simulations of this "limited imitation contagion" model on Poisson random 
graphs, finding agreement between the mean-field theory and stochastic simulations. 



I. INTRODUCTION 

Networks are an exploding area of research due to the 
recognition of their generality and ubiquity in physical, 
biological, technological, and social settings. Dynamical 
processes taking place on networks are now recognized as 
the most natural description for a number of phenomena. 
These include neuron behavior in the brain cellular 
genetic regulation Q, ecosystem dynamics and stability 
Q, and infectious diseases [1]. This last category, the 
study of biological contagion, is in many ways similar to 
social contagion, which refers to the spreading of ideas, 
fashions, or behaviors among people [a, 0] . This concept 
underlies the vastly important contemporary area of viral 
marketing, driven by the ease with which media can be 
shared and spread through social network websites. 

In this work, we present results for a very general 
model of networked map dynamics, motivated by models 
of social contagion. Each node has a "response function," 
a map which determines the state the node will take in 
response to the states of nodes in its neighborhood. Our 
model is a type of boolean network 0] , and it is closely 
related to models of percolation ^ and magnetism 
We focus on the derivation and analysis of dynamical 
master equations, both exact and mean- field approxima- 
tions, that describe the expected evolution of the system 
state. We also show how certain dense network limits 



lead to the convergence of the dynamics to the average 
response function map dynamics. 

We then apply our general theory to a particular lim- 
ited imitation contagion model [l^]. Nodes, representing 
people, act according to competing tendencies of imita- 
tion and non-conformity. One can argue that these two 
ingredients are essential to al l tren ds; indeed, Simmel, 
in his classic essay "Fashion" |l957l |. believed that these 
are the main forces behind the creation and destruction 
of fashions. Our model is not meant to be quantitative, 
except perhaps in carefully designed experiments, but it 
captures the features with which we are familiar: some 
trends take off and some do not, and some trends are sta- 
ble while others vary wildly through time. Our model is 
closely related to the seminal work of Schelling [l2| and 
Granovetter [l3l |. 

In Section HXl we define the general model and its de- 
terministic and stochastic variants. In Section IIIH we 
provide an analysis of the model when the underlying 
network is fixed. In Section llVl we develop a mean- field 
theory of the model on generalized random networks. In 
Section |Vl we consider the model on Poisson random 
networks with a specific kind of response function that 
reflects the limited imitation we expect in many social 
contagion processes. For this specific case, we compare 
the results of simulations and theory. Finally, in Sec- 
tion IVU we present conclusions and directions for further 
research. 
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II. GENERAL MODEL 

Let W = {V,E) be a network with = \V\ nodes, 
where V is the node set and E is the edge set. We let 
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A = A(^) denote the adjacency matrix; entry Aij is the 
number of edges from node j to node i. Assign each 
node i ^ V a response function fi : [0, 1] {0, 1}, and 
let x(0) G {0, 1}-'^ be the vector of initial node states. At 
time step t, each node i computes the fraction 



(1) 



of their neighbors in ^ who are active and takes the state 

x,{t + l) = f,{Mt)) (2) 

at the next time step. 

The above defines a deterministic dynamical system 
given a network and set of response functions. This is 
called a realization of the model [3] . Each node is in ei- 
ther the or 1 state; we refer to these as the off/inactive 
and on/active states, respectively. In the context of con- 
tagion, these would be the susceptible and infected states. 
With these binary states, our model is a particular kind 
of Boolean network ■ Note that each node reacts only 
to the fraction of its neighbors who are active, rather 
than the absolute number, and the identities of the in- 
put nodes do not matter. Each node's input varies from 
to 1 in steps of 1/ki, where ki = X^jLi is node z's 
degree (in-degree if is not a simple graph). 

In the rest of this Section, we will describe some vari- 
ations of the basic model which also differentiate our 
model from the Boolean networks extant in the litera- 
ture. This is mainly due to the response functions, but 
also the type of random network on which the dynamics 
take place, varying amounts of stochasticity introduced 
into the networks and response functions, and the possi- 
bility of asynchronous updates. 





Rewiring network 


Fixed network 


Probabilistic response 


P-R 


P-F 


Deterministic response 


D-R 


D-F 



TABLE I. The four different ways we implement the model, 
corresponding to differing amounts of quenched randomness. 
These are the combinations of fixed or rewired networks and 
probabilistic or deterministic response functions. In the ther- 
modynamic limit of the rewired versions, where the network 
and response functions change every time step, the mean-field 
theory (Sec. IIV|) is exact. 



are parameterized solely by two thresholds, (^on and c^off, 
so the distribution of response functions is determined 
by the joint density P{(j)on-,4>o{i)- Again, the specific net- 
work and response functions define a realization of the 
model. When these are fixed for all time, we have, in 
principle, full knowledge of the possible model dynamics. 
Given an initial condition x(0), the dynamics x(i) are 
deterministic and known for all t > 0. As for all finite 
Boolean networks 0] , the system dynamics are eventually 
periodic, since the state space {0, 1}^ is finite. 

With the introduction of noise, the system is no longer 
eventually periodic. Fluctuations at the node level allow 
a greater exploration of state space, and the behavior is 
comparable to that of the general class of discrete-time 
maps. Roughly speaking, the mean-field theory we de- 
velop in Section llVI becomes more accurate as we intro- 
duce more stochasticity. 

We introduce randomness in two parts of the model: 
the network and response functions. Allowing for the 
network and responses to be either fixed for all time or 
resampled each time step and taking all possible combi- 
nations yields four different designs (see Table [J) . 



A. The networks considered 

The mean-field analysis in Section HVl is applicable to 
any network which can be characterized by its degree 
distribution. The vast majority of the theory of ran- 
dom Boolean networks considers only regular random 
networks. Fortunately, such theories are easily general- 
ized to other types of networks with independently cho- 
sen edges, such as Poisson (Erdos-Renyi) and configura- 
tion model random networks 0, . We develop specific 
results for Poisson random networks, and these are the 
networks considered in Section IVl 



1. Rewired networks 

First, the network itself can change every time step. 
This is the rewiring (R), as opposed to fixed (F), net- 
work case. For example, we could draw a new network 
from Q{N, ki^vg/N) every time step. This amounts to 
rewiring the links while keeping the degree distribution 
fixed, and it is alternately known as a mean field, an- 
nealed, or random mixing variant as opposed to a fixed 
network or quenched model Q- 



Stochastic variants 



2. Probabilistic responses 



The specific network and response functions determine 
exactly which behaviors are possible. These are chosen 
from some distribution of networks, such as Q{N, fcavg/-^) 
(Poisson random networks on A'^ nodes with edge prob- 
ability ka.vg/N), and some distribution of response func- 
tions. In the example of Section|Vl the response functions 



Second, the response functions can change every time 
step. This is the probabilistic (P), as opposed to the de- 
terministic (D), response function case. For our social 
contagion example, there needs to be a well-defined dis- 
tribution P(0on, ^off) for the thresholds. This amounts to 
having a single response function, the expected response 
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function 

7(0) = / d0on / dc^off- -P(0on,0off)/((^;<^on,</>off)- (3) 



We call / : [0, 1] — ?► [0, 1] the probabilistic response func- 
tion. Its interpretation is the following. For an updating 
node with a fraction (p of active neighbors at the current 
time step, then, at the next time step, the node assumes 
the active state with probability f{(l)) and the inactive 
state with probability 1 — /{(/)). 



3. Temperature in the system 

In this paper, the network and response functions are 
either fixed for all time or resampled every time step. 
One could tune smoothly between the two extremes by 
introducing rates at which these reconfigurations occur. 
These rates are inversely related to quantities that behave 
like temperature (one for the network and another for the 
response functions). Holding the network or response 
functions fixed corresponds to zero temperature, since 
there are no fiuctuations. The stochastic and rewired 
cases correspond to high or infinite temperature, because 
reconfigurations occur every time step. 



C. Update synchronicity 

Finally, we introduce a parameter a for the probabil- 
ity that a given node updates. When a = 1, all nodes 
update each time step, and the update rule is said to be 
synchronous. When a ?a 1/N, only one node is expected 
to update with each time step, and the update rule is 
said to be effectively asynchronous. This is equivalent to 
a randomly ordered sequential update. For intermediate 
values, a is the expected fraction of nodes which update 
each time step. 



III. FIXED NETWORKS 

Consider the case where the response functions and 
network are fixed (D-F), but the update may be syn- 
chronous or asynchronous. Extend the definition of Xi{t) 
to now be the probability that node i is in the active state 
at time t. Note that this agrees with our previous defini- 
tion as the state of node i when Xi{t) = or 1. Then the 
Xi follow the master equation 



1^7 = 1 , 



x,{t + 1) = afi 



+ (4) 



which can be written in matrix-vector notation as 
x(t + 1) = af (rx(t)) + (1 - a)x(t). 



(5) 



Here T = D is sometimes called the transition prob- 
ability matrix (since it also occurs in the context of a 



random walker), D is the diagonal degree matrix, and 
f = ifz) 0- If a = 1, then x(t) G {0, 1}^ and we re- 
cover the fully deterministic response function dynamics 
given by Q and 



A. Asynchronous limit 

Here, we show that when a ?a time is effectively 
continuous and the dynamics can be described by an ordi- 
nary differential equation. This is similar to the analysis 
of Gleeson 16|. Consider Eqn. [SJ Subtracting x(t) from 
both sides and setting Ax(t) = x(t + 1) — x(t) and At = 1 
yields 



Ax(t) 
At 



a(f(rx(t))-x(t)). 



(6) 



Since a is assumed small, the right hand side is small, 
and thus Ax(t) is also small. Making the continuum ap- 
proximation dx(t)/dt f« Ax(t)/At yields the differential 
equation 



^=a(f(rx)-x) 



(7) 



The parameter a sets the time scale for the system. Be- 
low we see that, from their form, similar asynchronous, 
continuous time limits apply to the dynamical equations 
in the densely connected case, Eqn. ([8j), and in the mean- 
field theory, Eqns. ^ and (fT2]) . 



B. 



Dense network limit for Poisson random 
networks 



The following result is particular to Poisson random 
networks, but similar results are possible for other ran- 
dom networks with dense limits. The normalized Lapla- 
cian matrix is defined as £ = / — D~^/^ AD~^/'^ , where 
/ is the identity 0. So T = D-^/^{I - C)D^/^. By 
Oliveira [l^, when feavg is 0(log A^) there exists a typical 
Laplacian matrix = In — InI-Jj/N such that the 
actual C ~ in the induced 2-norm with high proba- 
bility [1 AT is the length- column vector of ones] . In this 
limit, if we approximate the degrees as uniform, ki ~ fcavg 
for ah i G V, then T « T^yp = l^lJ^/N. So T effectively 

averages the node states: Tx(t) J2iLi Xi{t)/N = 4>{t). 
Without a subscript, (f){t) denotes the active fraction of 
the network at time t. We make the above approxima- 
tion in Eqn. and average that equation over all nodes, 
finding 

(Pit + 1) = a/(0(t)) + (1 - a)0(t) = $(0(t); a), (8) 

where we have assumed that is large and the aver- 
age of nodes' individual response functions "^fLifi/N 
converges in a suitable sense to the stochastic response 
function /, Eqn. ([3]). This amounts to assuming a law 
of large numbers for the response functions, i.e., that the 
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sample average converges to the expected function. Note 
that a tunes between the probabilistic response function 
$(0; 1) = /((^) and the hue $(0;O) = (/). Also, the fixed 
points of $ are fixed points of /, but their stability will 
depend on a. When the network is dense, it ceases to 
affect the dynamics, since each node sees a large number 
of other nodes. Thus the network is effectively the com- 
plete network. In this way we recover the map models of 
Granovetter and Soo: 



3 way \ 
ng 



IV. MEAN-FIELD THEORY 

Making a mean-field calculation refers to replacing 
the complicated interactions among many particles by 
a single interaction with some effective external field. 
There are analogous techniques for understanding net- 
work dynamics. Instead of considering the interac- 
tions among the N nodes, network mean- field theories 
derive self-consistent expressions for the overall behavior 
of the network after averaging over large sets of nodes. 
These have been fruitful in the study of random Boolean 
network s |20| and can work well when networks are non- 
random [21|. 

We derive a mean-field theory, in the thermodynamic 
limit, for the dynamics of the general model by blocking 
nodes according to their degree class. This is equivalent 
to nodes retaining their degree but rewiring edges every 
time step. The model is then part of the well-known 
class of random mixing models with non-uniform con- 
tact rates. Probabilistic (P-R) and deterministic (D-R) 
response functions result in equivalent behavior for these 
random mixing models. The important state variables 
end up being the active density of stubs, i.e. half-edges 
or node-edge pairs. In an undirected network without 
degree-degree correlations, the state is described by a 
single variable p{t). In the presence of correlations we 
must introduce more variables {pkii)} to deal with the 
relevant degree classes. 



A. Undirected networks 

To derive the mean-field equations in the sim- 
plest case — undirected, uncorrelated random networks — 
consider a degree k node at time t. The probability that 
the node is in the active state at time t+ 1 given a density 
p of active stubs is 

Fk{p; f) = E f^)/''' (1 - P)'~' fij/k) ' (9) 

where each term in the sum counts the contributions from 
having 0, 1, . . . , fc active neighbors. Now, the probability 
of choosing a random stub which ends at a degree k node 
is qk = kpk I feavg in an uncorrelated random network . 
This is sometimes called the edge-degree distribution. So 



if all of the nodes update synchronously, the active den- 
sity of stubs at t + 1 will be 

oo oo , 



k=l 



k=l 



Finally, if each node only updates with probability a, we 
have the following map for the density of active stubs: 



p{t + l) = ag{p{ty,pk,f) + {l-a)pit) 
= G {p{t);pk,f, a) . 



ill) 



By a similar argument, the active density of nodes is 
given by 



cj){t + l) = ah{pit);pk,f) + il-a) (t){t) 
= H (p{t),(t){t);pkj,a) , 



where 



h {p; PkJ) = ^ PkFk {p; J) 

k=0 



(12) 



(13) 



Note that the edge-oriented state variable p contains all 
of the dynamically important information, rather than 
the node-oriented variable (h. 



B. Analysis of the map equation 

The function Fk{p; f) is known in polynomial approx- 
imation theory as the feth Bernstein polynomial (in the 
variable p) of / (2^ . Bernstein polynomials have im- 
portant applications in computer graphics due to their 
"shape-preserving properties" 23|. The Bernstein oper- 
ator Bfc takes f ^ Fk- This is a linear, positive operator 



which preserves convexity for all k and exactly interpo- 
lates the endpoints /(O) p and /(I). Immediate conse- 
quences include that each Fk is a smooth function and 
the A;th derivatives fI:^\x) — )■ f^^\x) where f^^\x) ex- 
ists. For / concave down, such as the tent or logistic 
maps, then Fk is concave down for all k and Fk increases 
to / {Fk /) as fe — > oo. This convergence is typically 
slow. Importantly, Fk J implies that g [p',PkT P) < / 
for any degree distribution pk ■ 

In some cases, the dynamics of the undirected mean- 
field theory given by p{t + 1) = G{p{t)), Eqn. (fTT]l . 
are effectively those of the map from the dense limit 
Eqn. ((8]). We see that g, Eqn. (fTOj) . can be seen as the 
expectation of a sequence of random functions Fk under 
the edge-degree distribution qk- Indeed, this is how it 
was derived. From the convergence of the -F^'s, we ex- 
pect that g{p;pk,f) ~ f{p) if the average degree fcavg 
is "large enough" and the edge-degree distribution has 
a "sharp enough" peak about fcavg (we will clarify this 
soon). Then as fcavg ~^ oo, the mean- field coincides with 
the dense network limit we found for Poisson random net- 
works, Eqn. [H Some thought leads to a sufficient condi- 
tion for this kind of convergence: the standard deviation 
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o'(^avg) of the degree distribution must be o(feavg)- In 
Appendix El we prove this as Lemma [T] 

In general, if the original degree distribution ph is char- 
acterized by having mean k^^g, variance a"^, and skewness 
7i, then the edge-degree distribution will have mean 
kavg + cr'^/ka.vg and variance cr^ [1 + 7i(7/A;avg - (cr/feavg)^]- 
Considering the behavior as fcavg — > oo, we can conclude 
that requiring a = o(fcavg) and 71 = o(l) are sufficient 
conditions on pk to apply Lemma [T] Poisson degree dis- 
tributions (cr = -^fcavg and 7i = fcavg^^) fit these criteria. 

C. Generalized random networks 

In more general random networks, nodes can have both 
undirected and directed incident edges. We denote node 
degree by a vector k = {k^^\ k^^\ k'^"^)^ (for undirected, 
in-, and out-degree) and write the degree distribution 
as = -P(k). There may also be correlations between 



node degrees. We encode correlations of this type by the 
conditional probabilities 

p^\.i = P(k, undirected|k') 
p^]^i = P(k, incoming|k') 
J'k°k/ = P(k, outgoing|k'), 

the probability that an edge starting at a degree k' node 
ends at a degree k node and is, respectively, undirected, 
incoming, or outgoing relative to the destination degree 
k node. We introduced this convention in a series of pa- 
pers [13, • These conditional probabilities can also be 
defined in terms of the joint distributions of node types 
connected by undirected and directed edges. We omit 
a detailed derivation, since it is similar to that in Sec- 
tion IIV Al and similar to the equations for the time evo- 
lution of a contagion process [2J, Eqns. (13-15)] [see also 



The result is a coupled system of equations for the density of active stubs which now may depend on node type (k) 
and edge type (undirected or directed): 
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The active fraction of nodes at a given time is 

<^(t + l) = (l-a)<^(t)+«E^"^E E 

ir n -n o .-n \ Ju J \ Ji J 
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(15) 



(16) 



Because these expressions are very similar to the undi- 
rected case, we expect similar convergence properties to 
those in Sec. IIVBI However, an explicit investigation of 
this convergence is beyond the scope of the current paper. 



V. LIMITED IMITATION CONTAGION MODEL 

As a motivating example of these networked map dy- 
namics, we study an extension of the classical threshold 
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models of social contagion [such asll2l. Il3l. l27H29l among 
others]. What differentiates our limited imitation conta- 
gion model from the standard models is that the response 
function includes an off-threshold, above which the node 
takes the inactive state. We assign each node i E V 
an on-threshold 4>on,i and an off-threshold 0off,ii requir- 
ing < (^on.i < 0off,i < 1- Node Vs response function 
M4>i) = M4>i',(pon,i',4>oS,i) is 1 if (/>on,i < 0i < 4>oS,i and 
otherwise. See Figure [T] for an example on-off threshold 
response function. 

This is exactly the model of Granovetter and Soong 
[T9| . but on a network. We motivate this choice with 
the following [also see[l^. (1) Imitation: the on state 
becomes favored as the fraction of active neighbors sur- 
passes the on-threshold (bandwagon effect). (2) Non- 
conformity: the on state is eventually less favorable with 
the fraction of active neighbors past the off-threshold (re- 
verse bandwagon, snob effect). (3) Simplicity: in the 
absence of any raw data of "actual" response functions, 
which are surely highly context-dependent and variable, 
we choose arguably the simplest deterministic functions 
which capture imitation and non-conformity. 

A crucial difference between our model and many re- 
lated threshold models is that, in those models, an acti- 
vated node can never reenter the susceptible state. Glee- 
son and Cahalane [2^ call this the permanently active 
property and elaborate on its importance to their anal- 
ysis. Such models must eventually reach a steady state. 
When the dynamics are deterministic, this will be some 
fixed fraction of active nodes. The introduction of the off- 
threshold builds in a mechanism for node deactivation. 
Because nodes can now recurrently transition between 
on and off states, the deterministic dynamics can exhibit 
a chaotic transient (as in random boolean networks 0]), 
and the long time behavior can be periodic with poten- 
tially high period. With stochasticity, the dynamics can 
be truly chaotic. 

The networks we consider are Poisson random net- 
works from Q{N, ka_yg/N). The thresholds c^on and (/)off 
are distributed uniformly on [0, 1/2) and [1/2, 1), respec- 
tively. This distribution results in the probabilistic re- 
sponse function (see Figure [2]) 
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FIG. 1. An example on-ofF threshold response function. 
Here, (pon = 0.33 and (f>oE = 0.85. The node will "activate" if 
(/'on < tp < f/)off, where is the fraction of its neighbors who 
are active. Otherwise it takes the "inactive" state. 
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FIG. 2. The tent map probabilistic response function /(p), 
Eqn. p7p . used in the limited imitation contagion model. 
This is compared to the edge maps g{p;ka,vg) = g{p;pk, f), 
Eqn. (|10p . shown for fcavg = 1, 10, 100. These pk are Poisson 
distributions with mean /cavg. As /cavg increases, q(^P] ^avg) 
increases to f{p)- 



A. Analysis of the dense limit 



ff < < 1/2, 
2cj) ff 1/2 < < 1. 



(17) 



The tent map is a well-known chaotic map of the unit in- 
terval [sOj. We thus expect the limited imitation model 
with this probabilistic response function to exhibit simi- 
larly interesting behavior. 



When the network is in the dense limit (Section IlII Bp . 
the dynamics follow (f){t + 1) = ^{(f){t); a), where 

^ f (1 + a)0 ff < < 1/2, (18) 

~ I (1 - 3a)(/» + 2a if 1/2 < < 1. 

Solving for the fixed points of $((^;a), we find one at 
(j) = and another at 4> = 2/3. When a < 2/3, the 
nonzero fixed point is attracting for all initial conditions 
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update probability a 

FIG. 3. Bifurcation diagram for the dense map $(0;q;), 
Eqn. (jlSp . This was generated by iterating the map at 1000 
a values between and 1. The iteration was carried out with 
3 random initial conditions for 10000 time steps each, dis- 
carding the first 1000. The (/>-axis contains 1000 bins and the 
invariant density, shown by the grayscale value, is normalized 
by the maximum for each a. With a < 2/3, all trajectories 
go to the fixed point at = 2/3. 



except (j) = 0. When a = 2/3, [1/2,5/6] is an interval 
of period-2 centers. Any orbit will eventually land on 
one of these period-2 orbits. When a > 2/3, this inter- 
val of period-2 centers ceases to exist, and more compli- 
cated behavior ensues. Figure [3] shows the bifurcation 
diagram for ^((/);a). From the bifurcation diagram, the 
orbit appears to cover dense subsets of the unit interval 
when a > 2/3. The bifurcation diagram appears like 
that of the tent map (not shown; see [lol. [30|) except the 
branches to the right of the first bifurcation point are 
separated here by the interval of period-2 centers. 



The effect of conformsits, an aside 

Suppose some fraction c of the population is made up 
of individuals without any off-threshold (alternatively, 
each of their off-thresholds (pos = 1). These individ- 
uals are conformist or "purely pro-social" in the sense 
that they are perfectly happy being part of the majority. 
For simplicity, assume a = 1. The map $(0; c) = 2(j) for 
0<(/)< l/2and2-2(l-c)(^forl/2<(/)< 1. Ifo 1/2, 
then the equilibrium at 2/3 is stable. Pure conformists, 
then, can have a stabilizing effect on the process. We 
expect a similar effect when the network is not dense. 



B. Mean- field 

Here we show how we compute the mean-field maps 
derived in Section IIVI In this specific example, we can 
write the degree-dependent map Fk{p; f) in terms of in- 
complete regularized beta functions Iz{a, b) Since / 
is understood to be the tent map, we will write (p; /) = 
Ffc(p). We find that 

Fkip) = 2p-4pIp{M,k~M), (19) 

where we have let M = [k/2\ for clarity ([-J and [•] 
are the floor and ceiling functions). The details of this 
derivation are given in Appendix [Bl 

The map g{p',Pk,f) is parametrized here by the net- 
work parameter fcavg; since pk is fixed as a Poisson dis- 
tribution with mean feavg and / is the tent map, and we 
write it as simply g{p;k^^g). To evaluate g{p;ki^vg), we 
compute Fk{p) using Eqn. (|19p and constrain the sum 
in Eqn. (fTUI) to values of k with [feavg — S-^fcavgJ < 
^ < [^avg + S-^feavgl • This computes contributions to 
within three standard deviations of the average degree 
in the network, requiring only 0(-\/fe^vg) evaluations of 
Eqn. (fT9|) . The representation in Eqn. (|19p allows for 
quick numerical evaluation of F^ip) for any k, which we 
performed in MATLAB. 

In Figure [21 we show g{p;kg^vg) for feavg = 1,10, 
and 100. We confirm the conclusions of Section IIV Bl 
5(p;^avg) is bounded above by f{p), and g{p]kscv^) / 
f{p) as fcavg -> oo. Convergence is slowest at p = 1/2, 
where the kink exhibited by the tent map has been 
smoothed out by the effect of the Bernstein operator. 

C. Simulations 

We performed stochastic simulations of the limited im- 
itation model for the D-F, P-F, and P-R designs, in the 
abbreviations of Table HI Unless otherwise noted, = 
10^. For all of the bifurcation diagrams, the first 3000 
time steps were considered transient and discarded, and 
the invariant density of p was calculated from the follow- 
ing 1000 points. For plotting purposes, the invariant den- 
sity was normalized by its maximum at those parameters. 
For example, in Figure [3] we plot P((/!)|a)/ max^ P((/)|a) 
rather than the raw density P{(j)\a). 

To compare the mean-field theory to those simula- 
tions, we numerically iterated the edge map p{t + 1) = 
G {p{t); kgjvg, a) for different values of a and fcavg- We 
then created bifurcation diagrams of the possible behav- 
ior in the mean- field as was done for the simulations. 



D. Results 

To provide a feel for the deterministic dynamics, we 
show the result of running the D-F model on a small net- 
work in FigureHl Here, N = 100 and fcavg = 17. Starting 
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time 



FIG. 4. Deterministic (D-F) dynamics on a small network. 
Here, N = 100 and fcavg = 17. On the left, we plot the 
state evolution over time. The upper plot shows individual 
node states (black = active) sorted by their eventual level of 
activity, and the lower plot shows the total number of active 
nodes. We see that the contagion takes off, followed by a 
transient period of unstable behavior until time step 80, when 
the system enters a macroperiod-4 orbit. Note that individual 
nodes exhibit different microperiods (see Sec. IVD|I . On the 
right, we show the network itself with the initial seed node in 
black in the lower right. 



from a single initially active node at t = 0, the active pop- 
ulation grows monotonically over the next 6 time steps. 
From t = 6 to t = 80, the transient time, the active frac- 
tion varies in a similar manner to the dynamics in the 
stochastic and mean- field cases. After the transient, the 
state collapses into a period-4 orbit. We call the overall 
period of the system its "macroperiod," while individual 
nodes may exhibit different "microperiods." Note that 
the macroperiod is the lowest common multiple of the 
individual nodes' microperiods. In Figure [H we observe 
microperiods 1, 2, and 4 in the timeseries of individual 
node activity. In other networks, we have observed up 
to macroperiod 240 [l^]- A majority of the nodes end 
up frozen in the on or off state, with approximately 20% 
of the nodes exhibiting cyclical behavior after collapse. 
The focus of this paper has been the analysis of the on- 
off threshold model, and the D-F case has not been as 
amenable to analysis as the stochastic cases. We offer a 
deeper examination through simulation of the determin- 
istic case in • 

We explore the mean-field dynamics by examining the 
limiting behavior of the active edge fraction p under the 
map G (p; fcavg, a). We simulated the map dynamics for 
a mesh of points in the (fcavg,") plane. We plot the 3- 
dimensional (3-d) bifurcation structure of the mean-field 
theory in FigureO We also show 2-d bifurcation plots for 
fixed fcavg and a slices through this volume in Figures [6] 
and [71 For more visualizations of this bifurcation struc- 
ture, see Appendix O In all cases, the invariant density 
of p is normalized by its maximum for that (fcavg, ") pair 
and indicated by the grayscale value. 




FIG. 5. The 3-dimensional bifurcation diagram computed 
from the mean-field theory. The axes X = average degree fcavg, 
Y = update probability a, and Z = active edge fraction p. The 
discontinuities of the surface are due to the limited resolution 
of our simulations. See Figure[6]for the parameters used. This 
was visualized using Paraview, and a file is available in the 
online Appendix. 



The mean-field map dynamics exhibit period-doubling 
bifurcations in both parameters fcavg and a. Visualiz- 
ing the bifurcation structure in 3-d (Figure E]) shows in- 
terlacing period-doubling cascades in the two parame- 
ter dimensions. These bifurcations are more clearly re- 
solved when we take slices of the volume for fixed pa- 
rameter values. The mean- field theory (Figure [6]) closely 
matches the P-R simulations (Figure[7]). The first deriva- 
tive ^(p; fcavg, a) < |if for any finite fcavg, so the 
bifurcation point a = 2/3 which we found for the dense 
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average degree k^^^ 

FIG. 6. Mean-field theory bifurcation diagram slices for various fixed values of fcavg and a. The top row (a-c) shows slices for 
fixed fcavg- As fcavg — >■ c», the /cavg-slice bifurcation diagram asymptotically approaches the bifurcation diagram for the dense 
map, FigurelHl Note that the first bifurcation point, near 2/3, grows steeper with The bottom row (d-f) shows 

slices for fixed a. The resolution of the simulations was a = 0.664, 0.665, . . . , 1, k^vg = 1, 1.33, . . . , 100, and p bins were made 
for 1000 points between and 1. 




(b) k =50 



0.8 0,9 1 




0.8 0.9 1 

update probability a 




(d) a = 0.68 



(e) a = 0.92 



(f) a = 1 
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average degree k 
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avg 



FIG. 7. Bifurcation diagram from fully stochastic (P-R) simulations, made in the same way as Figure [S] The bifurcation 
structure of these stochastic simulations matches that of the mean- field theory (Figure [6]), albeit with some blurring. 
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map $ is an upper bound for the first bifurcation point 
of G. The actual location of the first bifurcation point 
depends on fcavg, but a = 2/3 becomes more accurate 
for higher feavg (it is an excellent approximation in Fig- 
uresEt and[7t, where fcavg = 100). When a = 1, the first 
bifurcation point occurs at feavg ^ 7. 

The bifurcation diagram slices resemble each other and 
evidently fall into the same universality class as the lo- 
gistic map [HjIsI]. This class contains all 1-d maps with 
a single, locally-quadratic maximum. Due to the prop- 
erties of the Bernstein polynomials, Fh{p; f) will univer- 
sally have such a quadratic maximum for any concave 
down, continuous / [22|. So this will also be true for 
5(/o;^avg;/) with fcavg finite, and we see that fcavg par- 
tially determines the amplitude of that maximum in Fig- 
ure [21 Thus fcavg acts as a bifurcation parameter. The 
parameter a tunes between G (/o; fcavg, 1) = 9{p', k^^g, f) 
and G {p;kavg,0) = p, so it has a similar effect. Note 
that the tent map / and the dense limit map ^ are 
kinked at their maxima, so their bifurcation diagrams 
are qualitatively different from those of the mean-field. 
The network, by constraining the interactions among the 
population, causes the mean-field behavior to fall into a 
different universality class than the individual response 
function map. 



VI. CONCLUSIONS 

We have described a very general class of synchronous 
or asynchronous, binary state dynamics occuring on net- 
works. We obtained an exact master equation and 
showed that, when random networks are sufficiently 
dense, the networked dynamics approach those of the 
fully-connected case. We developed a mean-field theory 
and found that it also predicted the same limiting behav- 
ior. The convergence of the mean-field map to the aver- 
age response function is related to the Bernstein poly- 
nomials, allowing us to employ many previous results in 
order to analyze the mean-field map equation. We also 
extended those mean-field equations to correlated ran- 
dom networks. 

The general model we describe was motivated by the 
the limited imitation model of social contagion. We see 
that including an aversion to total conformity results in 
more complicated, even chaotic dynamics, as opposed to 
the simple spreading behavior typically seen in the single 
threshold case. The theory developed for the general case 
successfully captured the behavior of the stochastic net- 
work dynamics. We have focused on the rich structure of 
bifurcations as the two parameters, update synchronicity 
a and average degree feavg, were varied. We see that the 
universality class of the dynamics matches those of the 
logistic map. Using the mean-field theory, we can un- 
derstand this as a result of the smoothing effect of the 
Bernstein polynomials on the tent map average response 
function. However, this universality class will appear for 
any unimodular, concave down response function. 



The deterministic case, which we have barely touched 
on here, merits further study [see [l3|. In particular, 
we would like to characterize the distribution of periodic 
sinks, how the collapse time scales with system size, and 
how similar the transient dynamics are to the mean-field 
dynamics. 

Furthermore, the model should be tested on realistic 
networks. These could include power law or small world 
random networks, or real social networks gleaned from 
data. In a manner similar to Melnik et al. [2l|, one could 
evaluate the accuracy of the mean-field theory for real 
networks. 

Finally, the ultimate validation of this model would 
emerge from a better understanding of social dynamics 
themselves. Characterization of people's "real" response 
functions is therefore critical [some work has gone in this 
direction; see H, H, [H, [S^l . Comparison of model output 
to large data sets, such as observational data from social 
media or online experiments, is an area for further experi- 
mentation. This might lead to more complicated context- 
and history-dependent models. As we collect more data 
and refine experiments, the eventual goal of quantifiably 
predicting social behavior, including fashions and trends, 
seems achievable. 
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Appendix A: Proof of Lemma 1 

Lemma 1. For k > 1, let fk be continuous real-valued 
functions on a compact domain X with fk ^ f uni- 
formly. Let pk be a probability mass function on Z"*" 
parametrized by its mean p and with standard deviation 
(t{p), assumed to be o{p). Then, 

lim y^Pkfk = /. 

\k=0 / 

Proof. Suppose < a < 1 and let K = [p — p°-\ . Then, 

CXD K oo 

9 = ^Pkfk = ^Pkfk+ ^ Pkfk- (Al) 

fc=0 fc=0 k=K+l 

Since fk^f uniformly as — > oo, for any e > we can 
choose p large enough that 

\fk{x)-f{x)\<e (A2) 
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for all k > K and all x 6 X. Without loss of generality, 
assume that \fk\ < 1 for all k. Then, 



|5-/I< 



The o"/^" term is a consequence of the Chebyshev in- 
equality [13] applied to the first sum in (|A1|) . Since 
a grows sublinearly in /i, this term vanishes for some 
< a < 1 when we take the limit fi ^ oo. The e term 
comes from using (|A2p in the second sum in ()Aip . and it 
can be made arbitrarily small. □ 



form of Eqn. (fT71) to write 

M 



ft(P)^EG) 

k 

+ E 



p^"(l-p)'^-^" 2 



El 

A: 



2-2p-2^r p^(l-^)'=-^- 



We have used the fact that the binomial distribution 
— p)^~^ sums to one and has mean kp. For 
n < AI, we have the identity 

AI 



3=0 



p^il-pf-^ = p^{k)^h_p{k-M,M-n + l) 

(B2) 

where Ix{ci, b) is the regularized incomplete beta function 
and (fc)n = k{k~ 1) • • • (fc — (n— 1)) is the falling factorial 
[Sll . ISq . This is an expression for the partial (up to M) 
nth factorial moment of the binomial distribution with 
parameters k and p. Note that when n = we recover 
the well-known expression for the binomial cumulative 
distribution function. From Eqns. (|B1|) and (|B2[) . we 
arrive at Eqn. (fT9]) . 



Appendix C: Online material 



Appendix B: Beta function representation of Fk 



We now show how, when / is the tent map (jT7|) . the 
degree-fc map F^^p; f) can be written in terms of incom- 
plete regularized beta functions. First, use the piecewise 



To better explore the 3-d mean field bifurcation 
structure, we created movies of the the fcavg and a 
slices as the parameters are dialed. We also provide 
a VTK file with the 3-d bifurcation data, view- 
able in Paraview or other software. Videos of the 
individual-node dynamics for small networks in the 
D-F and P-F cases are shown for some parameters 
which produce interesting behavior. The D-F, P-F, 
and P-R cases were implemented in Python, and we 
provide our code. All of the above is available from 
,http : / / www . uvm . edu/ storylab/ share/ paper s/harris2013a/[ 



[1] S. Coombes, B. Doiron, K. Josic, and E. Shea-Brown, 
Philosophical Transactions of the Royal Society A: Math- 
ematical, Physical and Engineering Sciences 364, 3301 
(2006). 

[2] R. Milo, S. Shen-Orr, S. Itzkovitz, N. Kashtan, 

D. Chklovskii, and U. Alon, Science 298, 824 (2002). 
[3] J. Bascompte, Science 325, 416 (2009). 



[4] P. Rohani, X. Zhong, and A. King, Science 330, 982 
(2010). 

[5] J. Ugander, L. Backstrom, C. Marlow, and J. Kleinberg, 
PNAS 109, 5962 (2012). 

[6] D. Centola, Science 329, 1194 (2010). 

[7] M. Aldana, S. Coppersmith, and L. P. KadanofF, in 
Perspectives and Problems in Nonlinear Science, edited 
by E. Kaplan, J. E. Marsden, and K. R. Sreeni- 



12 



vasan (Springer, New York, 2003), chap. 2, pp. 23-90, 
[http : //arxiv . org/abs/nl in/0204062 , 
[8] D. Stauffer and A. Aharony, Introduction to Percolation 

Theory (Taylor and Francis, London, 1994), 2nd ed. 
[9] M. E. J. Newman, SIAM Review 45, 167 (2003). 
[10] P. S. Dodds, K. D. Harris, and C. M. Danforth (2013), 

http : //arxiv . org/abs/1208 . 0255 
[11] G. Simmel, American Journal of Sociology 62, 541 
(1957). 

[12] T. C. Schelling, Journal of Mathematical Sociology 1, 
143 (1971). 

[13] M. Granovetter, American Journal of Sociology 83, 1420 
(1978). 

[14] B. BoUobas, Random Graphs (Cambridge University 

Press, 2001), 2nd ed. 
[15] When there are isolated nodes, those nodes' degrees are 

0, and thus the denominator in the master equation is 

zero and undefined. If the initial network contains 

isolated nodes, we set all entries in the corresponding 

rows of T to zero. 
[16] J. P. Gleeson, Phys. Rev. E 77, 046117 (2008). 
[17] D. B. West, Introduction to Graph Theory (Prentice Hall, 

Upper Saddle River, NJ, 2001), 2nd ed. 

[18] R. I. Qliveira (2010), [http://arxiv.org/abs/0911.0600l 

0911.0600. 

[19] M. Granovetter and R. Soong, Journal of Economic Be- 
havior and Organization 7, 83 (1986). 

[20] B. Derrida and Y. Pomeau, Europhys. Lett. 1, 45 (1986). 

[21] S. Melnik, A. Hackett, M. A. Porter, P. J. Mucha, and 
J. P. Gleeson, Phys. Rev. E 83, 036112 (2011). 

[22] G. M. Phillips, Interpolation and Approximation by Poly- 
nomials (Springer, 2003), see chapter 7. 



[23 
[24 
[25; 

[26; 



J. M. Peha, ed., Shape Preserving Representations in 
Computer- Aided Geometric Design (Nova Science Pub- 
lishers, Inc, 1999). 

J. L. Payne, K. D. Harris, and P. S. Dodds, 
Phys. Rev. E 84, 016110 (2011), URL 
http : //link . aps . org/doi/10 . 1103/PhysRevE . 84 .016110] 
P. S. Dodds, K. D. Harris, and J. L. Payne, 
Phys. Re v. E 83, 056122 (2011), URL 
|http : //link . aps . org/doi/10 . 1 103/PhysRevE . 83 . 056122( 



J. P. Gleeson and D. J. Cahalane, 
Phys. Rev. E 75, 056103 (2007), URL 
|http : //link . aps . org/doi/ 10 . 11 03/PhysRevE . 75 . 056103| 

[27] T. C. Schelling, Journal of Conflict Resolution 17, 381 
(1973). 

[28] D. J. Watts, PNAS 99, 5766 (2002). 
[29] P. S. Dodds and D. J. Watts, Phys. Rev. Lett. 92, 218701 
(2004). 

[30] K. T. AUigood, T. D. Sauer, and J. A. Yorke, Chaos: An 
Introduction to Dynamical Systems (Springer, 1996). 

[31] NIST, Digital Library of Mathematical Functions (2012), 
http : // dlmf . nist . gov 

[32] M. J. Feigenbaum, Journal of Statistical Physics 19, 25 

(1978) . 

[33] M. J. Feigenbaum, Journal of Statistical Physics 21, 669 

(1979) . 

[34] D. Centola, Science 334, 1269 (2011). 

[35] D. M. Romero, B. Meeder, and J. Kleinberg, in Proc. 20th 
ACM International World Wide Web Conference (2011). 

[36] R. L. Winkler, G. M. Roodman, and R. R. Britney, Man- 
agement Science 19, 290 (1972). 



